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ABS Ru, 


In airborne Anti-Submarine Warfare operations there is 
a Liteon requirement for maintaining an accurate relative 
plot of the sonobuoys with respect to the aircraft. This 
Study proposed a method for locating sonobuoys in a pat- 
ern Using aircraft-to-buoy slant range information. The 
method did not use trtangulation procedures and attempted 
to minimize the restrictions placed on the aircraft. The 


Study showed the feasibility of the proposed methodology and 


the approximate errors to be encountered, 
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E INTRODUCTION 


The fixed wing dincraft plays ean inpog santero le Simamecnlies 
borne antisubmarine warfare ms metod ay mine mips iiicl Gee iGsteiin 
cal employment of the aircraft is through ihe use of active 
and passive electronic sensors, the most useful of which is 
the sonobuoy. Sonobuoys are placed by the atrerat t. tnavarer- 
ous geometric patterns in the water as the tactical situa- 
tion predicts. The primary phases in the prosecution of a 
menitact are search, detection, classification, localization 
and attack. or in the rare case of a visual sighting of 
a submarine, the first four phases are almost exclusively 
conducted by the use of passive and active sonobuoys which 
transmit data, analyzed by sensor equipment on board the 
aircraft. This information is then synthesized by the air- 
craft tactical commander to yield an assessment of the tac- 
tical situation. 

One major problem in airborne ASW is the lack of accu- 
rate aircraft/sonobuoy referencing. Surface wind, sea cur- 
rents and errors in navigation quickly disperse the actual 
buoy pattern location from its assumed location. Moreover 
the buoys within the pattern frequently drift in a random 
manner, particularly anen the distance between buoys is 
great. lt its necessary to periodically update the actual 
buoy positions by a visual or electronic "mark on top.' 
ine eee entity of the aircraft is obviously reduced by this 


requirement.- In the fine localization phase where a combi- 


nation of two or more buoys in close proximity to the target 








is providing the a@ireraft witiiet ) data cece sae, 

to fly outbound to the tix “after spassing Oven munemG loses: 
buoy providing data. This tactic is necessary to provide 
the aircraft with final, accurate buoy location information 
must prior to the attack.  Frequentemissed i iiswan ec ot yi 
buted to attacking a fix position generated by buoys whose 
meeurate locations with respect to the aircratt were not 
known or recently updated. 

A better method for determining the location of sono- 
buoys is required. A method for locating buoy coordinates 
relative to the aircraft, remote from the buoy, would elimi- 
nate the need for ''on top" repositioning. Such a method 
would also offer the additional aircraft flexibility re- 
quired for the eventual increase in the tactical load placed 
on the crew. 

A system under development [Ref. 1,2] will provide for 
a slant range measurement to the buoy. This slant range 
would be more accurate than the aircraft navigation system 
and could be converted to horizontal range knowing the air- 
craft altitude. The range to a buoy can be taken at any 
time after it has been placed in the water up to the life 
of the buoy and in ste tact real environment providing the 
aircraft fs at sufficient altitude to receive the signal. 
The use of the ranges provided by the buoys under the pro- 
posed system [Ref. 1] serves as the subject of this study. 


While triangulation methods obviously provide a solu- 


Pion to the buoy location problem {Fig. 1] th alreratur 
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A/C Track 


horizontal range from A/C to buoy 


A/C location on track when first range 
measured 


A/C location on track when second range 
measured 


— 


separation angle between range measurements 


FIGURE I 


STANDARD TRIANGULATION PROCEDURE 
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track, line AB, between the buoy rangings required to com- 
plete triangulation should be approximately equal to the 
range of the buoy. The buoy fix position is determined by 
the intersection of the range arcs to the bu OY from the 
aircraft position. The errors realized in the buoy fix 
Bosition are, for the most part, due to inaccuracies tn the 
ranging system and inaccuracies in navigation from points 
A to B. If the leg AB [Fig. 1] representing the aircraft 
track is allowed to be in the order of length equal to the 
mime to the buoy such that the range arcs itntersect as 
close to perpendicular as possible, then the errors due to 
navigation and range measurement would be transferred to the 
buoy fix error in a nearly one-to-one ratio. If, however, 
the aircraft track, leg AB, between range measurements is 
short, inferring an angle 6 of only a few degrees, then the 
range arcs will intersect nearly parallel. The bearing 
error of the resulting buoy fix position may be extremely 
large. - 
tn that the length of the aircraft track required for 
optimal error reduction in triangulation necessarily re- 
Stricts the flexibility of the aircraft, this study wil 
focus on a method of inne the error by statistical 
methods using a rapid series of rangings to the buoy over 
a short interval. It is felt that the statistical methods 


of error reduction will allow an improved solution to the 


problem over standard triangulation. 











[b. (NATGORE ORS Gear Rosen 


The major portion of this ctudyeis centered atouma 
(1) determining a method for locating sonobuoys relative to 
mie aircraft, (2) feastbility of the systemeyith respect to 
the deployment of the aircraft in a tactical environment, 
and (3) the errors anticipated and how sensitive the pro- 
Seawire is to these errors. A residual effort conducted in 
conjunction with the fev anleenent of the sonobuoy location 
maomlen provides qa conceptual solution method for deter- 
minting an accurate aircraft ground speed and true course. 
iivesdevelopment concepts for the ground speed and true 
course determination are not complicated and should therefore 
be adapted easily to any aircraft in the fleet today which 
would carry the associated sonobuoy ranging equipment. 

Range information from the buoys will be available con- 
tinuously, sampled as destred. The slant range will be 
Eonverted to horizontal range as described by Ellis [Ref. l, 
Fig. 4]. There is no inherent ambiguity problem such as 
mirror images of fix positions and the solution is not 
quadrant or half plane restricted. Minor restrictions are 
placed on the aircraft for short periods of time. The air- 
Craft, in most cases, will be required to fly Straight and 
level at a con:tant true airspeed while sampling ranges to 
the buoys; however these intervals are quite short, four to 
menty seconds as the tactical situation dictates. This 


Bestriction is not considered to be excessive as wil] be 


explained tn the formal development of techniques for 
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sonobuoy location and ground speed determination. gene ets 
tain cases the measured ranges to the sonobuoys will be the 
only input required by the ye ecneaien generates a very 
appealing aspect to the solution technique. “In other cases 
the aircraft eoorclneece as well as the measured range in- 
formation will be required as the only inputs. 

It should be noted here that navigational errors pro- 
ducing erroneous geographical track of the aircraft are not 
critical in the problem solution as such. The primary pro- 
blem being attacked by this study is that of locating the 
sonobuoy relative to the aircraft. If, in fact, the DR 
track is erroneous then the positions determined for the 
buoys will be erroneous; however, their positions relative 
to the aircraft will be accurate as they will be plotted 
relative to the aircraft DR position. Accurate geographical 
track and sonobuoy location information can be updated sub- 
Sequent to the attack phase of the problem. Experience has 
Shown that two aircraft in the process of relieving each 
other may have a relative error of up to 15 miles in geo" 
Bepnical DR tracks. Therefore the transfer of geographical 
Peertions of buoys within a pattern_is frequently unsatis- 
factory. When it is necessary to turn the problem over to 
a relieving aircraft, tt will be easier to transfer a rela- 
tive picture, that is, range and bearing of all other buoys 
in the pattern from one reference buoy Sioee teowtine Teli ey img 
Empccraft. 


Although the P3-C aircraft will far surpass--any existing 
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ASW airborne weapon system in the fleet today it is not 
anticipated that the models used currently wil! be phased 
out in the near future. lt is important to consider devel- 
oping a system which may, with minimal modification, be 
adapted to sineae models of ASW aircraft. The ground speed 
and true course determination is developed with this in mind. 
Considering the contribution of locating sonobuoys from re- 
mote positions to the prosecution of a contact, it may be 
well worth while to consider the installation of a limited 


purpose computer in the current model VP aircraft vice more 


costly alternatives. 


il 





Lil. BUOY RANGING CONCEPT 


The sonobuoy has a transmitter built in for the purpose 
of transmitting to the aircraft the sensor signals received 
from the sonobuoy hydrophones. The carrier frequency asso- 
ciated with this signal is channelized in each buoy for 
buoy discrimination. A buoy ranging signal described tn 
detail by Ellis [Ref. 1] will be transmitted over the exist- 
ing aircraft receivers. No new receivers or transmitters 
for the signal will be required. A modification to the 


ey 


detector section of Ave cocoon receiver will be required 
for signal separation. 

The sonobuoy ranging signal is generated by a crystal 
oscillator placed in the buoy. The method of range deter- 
mination ts by comparison of the phase of the buoy signal 
with that of a standard frequency generated on board the 
atrcraft. The phase fae coencs is Converted “toms lant. tanece 

The sonobuoy undergoes a warmup period in the aircraft 
prior to drop. tt generates a stable frequency, ree and is 
compard with the aircraft frequency, ee The phase dif- 
ference at this time and at zero distance is measured, AG. 
At any time, t, after the sonobuoy is dropped in the water 
a new comparison is eee ns and fr and the phase 


AG, is determined. The resulting phase difference, AO, 


where AO. = AG. = AG. is then converted to slant range; 


~e 


Subsequently to horizontal range. The buoy range frequency 
has a wavelength of approximately 75 kilometers. One cycle 


difference in phase of the buoy signal from the reference 
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is equivalent to 75 km. distance Seis buey eae cm cigiemon 

10 km. and 85 km. would yield the same phase difference 
causing an ambiguity. lt must be known what multiple of 

75 km. the buoy is from the aircraft. It is not anticipated 
that this fr pcan will cause erroneous ranging. 

The range errors generated by the system are covered 
iedetail by Ellis [Ref. 1] but one error of primary con- 
cern in the tactical employment of the system is the linear 
increase in frequency with respect to time within the buoy 
oscillator. Gampilene frequency stabilization is not yet a 
"state-of-the-art'' realization. Testing to date with the 
oscillator currently in use in the buoy has generated an 
average drift rate comparable to a range error of 1040 meters 
per hour. tf a means was devised for the subtraction of the 
linear drift rate as a function of time, the remaining 


~ 


errors in the system could be aggregated to a classification 
of normal with a mean of zero and a variance of o.. Tas 
aggregation of errors being N(O,0,) would allow a broad 
application of statistical methods for the solution of the 
error reduction problem in the tactical employment of the 
buoy range. _It is therefore assumed that in the following 
study the average ies ocaron rate is eliminated from the 


problem and the assumption that the remaining errors are 


2 
N(0,0.) will motivate the development. 








IV. SUMMARY OF PROBEEM SSO EUT OnE Caio 


The following methods are developed for the solution 

of (1) Sonobuoy Location, (2) Aircraft Ground Speed and (3) 
Aircraft True Course. This section 5 ON eae non-techni = 
cal summary of the problem geometry, basic solution approach 
and the assumptions and restrictions governing the technique 
used in each of the three problem areas mentioned above. It 
is designed as a background preface to the more detailed 
mathematical development and should by itself acquaint the 
reader with a complete broad picture of what has been accom- 


plished. The mathematical development will be included in 


a subsequent section. 


we SONOBUOY LOCATION 

lt is desired to provide a system whereby the aircraft 
may, in any phase of the tactical prosecution of a contact, 
locate any or all of the sonobuoys in the pattern with re- 
spect to the aircraft. To provide the aircraft with the 
flexibility needed in the more acute phases of the problem, 
the system must not place heavy restrictions on the maneuver- 
ability of the aircraft such as is required in standard 
meiangulation techniques. A mere o yi teatenn the aircraft 
could obtain locations from a single position and time or 
with a minimum advance leg would certainly be preferred. 

Recalling the discussion in the last paragraph of the 
iatroduction, an extended aircraft track, line AB [Fig. 1], 


would have to be flown in order that 98 be sufficiently large 


to reduce fix error of the buoy position. In a tactical 
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situation it may not be feasible for the aircraft to fly 
the leg AB. lf the aircraft desired to maneuver between 
points A and B it must do so at the poss Dllecostmen mcd cdi 
tional navigational errors as aircraft navigation becomes 
degraded in many cases as a result of maneuvering. 

The geometry in Figure 2 was considered as a starting 
point for the development for sonobuoy location. This basic 
meometry will prevail in the description of the problem. 
Consider an aircraft with a buoy in the relative position 


own. The aircraft eceamhired to hold any true heading 
and airspeed for as long as it its required to measure 

ranges Ry, R2 and R3 [Fig. 3]. The time interval between 
each range measurement is short, approximately five seconds. 
Pee2z00 kts. an aircraft travels 111 yds./second. As many 
ranges as desired may be measured but more time is required 
and the additional secure obtained by more ranging tis not 
significant as a trade off in a critical tactical situation. 
This will be a tactical decision. The study uses three 
mage measurements throughout the descriptive phase. Addi- 
Mei ranges were used in testing. [It its assumed that in 
the short period of time the aircraft is required to fly 
straight and at a eeneeontt air speed that the coordinates 
(A, ,B.) [Fig. 3] are accurate relative to each other. Any 
mmaccuracies occuring wil] be absorbed in final fix error. 
The men aes R. and the aircraft coordinates (A, ,B.) are in- 
puts and are assumed to be known. The buoy coordinates 


(X,Y) are unknown. 


We 








CPA 





A/C Heading 


R. = horizontal range from A/C to buoy 


CPA = closest point of approach to buoy 
with respect to the A/C 


— 


6 = relative bearing of buoy with respect 
to the A/C 
FIGURE 92 


BASIC TWO DIMENSIONAL PROBLEM GEOMETRY 
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(A, .B,) 


(ey ) 


aaw oes he Se os O88 Se ee 


Gite Bia) H 


horizontal 


A/C coordinates at time of measurement of 
bel eee 


range R, to buoy 


coordinated of buoy (unknown) 


FIGURE 3 


SONOBUOY LOCATION PROBLEM GEOMETRY 
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a ee 


eading 


(xe) 


range to buoy as measured by A/C 
io ene 
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The model used to determine the (X,Y) position of the 
buoy is the method of Least Squares. The basic assumption 
fori the use of this model is that the errors are N(0,0.). 
mine approach is to minimize the sum ofthe squares of the 
errors between the measured ranges to the buoy and the 


actual ranges. The equation for the Sum of Squares; 


a > (0.-E.) where; 0.=observed value 
E.=expected value 


0.-E.=e, (error term) 


aes) 


i=1 


is adapted to the sonobuoy ranging problem; 


sf =) Coane where; R.=measured range 
Ri =expected value 


=[(X-A,) + (¥-B,)°]? 


The approach of minimizing the sum of the error squared 
terms “is that of evaluating the first derivative of SS at 


zero with respect to the X and Y values desired; 


n eC 2 
ss = J (e.)> = J {R, - £(x-a.)° + (¥-B,) 94} 
it i=] 
BSE OSS 7 
mie «Sy ee 


This development produces a set of normal equations from 


which the values of X and Y are determined by an iterative 
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process. 

The iterative process requires that an estimate of the 
buoy angle relative to the aircraft, 90, be supplied. This 
initial angle, 69, provides a starting point for ERemp ir o> 
cess which by making alternately successive approximations 
of X, Y and 6 the process converges on the values of X and 
Y with the minimum error possible for the geometry. Values 
Sie oo Up to and inctuding + 15° in error of the actual 
angle 0 were tried [Fig. 4]. There was little effect on 
the rapid eonseGenee 5 X and Y to within acceptable values. 
In consideration of the latitude in choosing a value of 69 
it is felt that the bearing supplied to the buoy by the DME 
indicator is satisfactory as an input value. The value of 
69 was considered to have been chosen as an estimate of the 
8 for the mic’-range measurement. 

Two studies were conducted to determine the feasibility 
of the Least Squares model for the solution of the sonobuoy 
location problem. The first study was a sensitivity analysis 
to determine the effect of the following on the convergence 
rate of the process; (1) distance to the buoy, (2) relative 
mmole of the buoy, 6, (3) values of 89 as an estimate of 
and (4) number of Pema cae in the solution of the pro- 
blem. 

The second study was an error analysis to examine the 
Beaect of range error on the final values of buoy coordi -— 


nates, X and Y. Errors of varying magnitudes were applied 


to the ranges in different combinations at different relative 
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Zo Buoy £} 
ae . 
7 +150 
J 
2 
| A/C 
A Heading 
A = A/C location when horizontal range to buoy 
measured. - 
6 = actual relative angle of buoy with respect 
to A/C heading 
6 = estimated values of relative angle of buoy; 


error applied in 1° increments up to + 15° 


: of actual relative angle 


FIGURE 4 


RELATIONSHIP BETWEEN ESTIMATED AND ACTUAL VALUES OF RELATIVE 
BUOY ANGLE WITH RESPECT TO A/c HEADING 


os 
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buoy angles. The results of these studies are explained in 


detail in Appendix A. 


B. GROUND SPEED DETERMINATION , 

The development of the method of ground speed deter- 
mination is based on the concept that if an aircraft was 
approaching a buoy (0°relative bearing) which was providing 
accurate range information then the change in the horizontal 
range to the buoy with respect to time would be, exactly, 
the ground speed. This is a fast, accurate and simple way 
to determine ground speed if it could be ascertained that 
the relative bearing was exactly 9° relative. This is 
difficult to do as the DME indicator does not provide that 
accuracy. As the relative angle of the buoy increases from 
zero the first derivative of the range to the buoy with time 
becomes a function of the cosine of the relative angle. If 
there existed a scheme for determining the change in range 
to the buoy with respect to time with the relationship of 
the relative angle then the qureundespeed  couldsbe determ tuned 
for any relative angle for the buoy. 

In Figure 5 a geometry similar to that of the previous 
section is used. Consider an aircraft approaching a buoy 
as indicated. Were the aircraft to continue on this head- 
ing it would pass the buoy at its closest point of approach 
feaen) at a point perpendicular to the aireraft track. “In 
the problem solution the aircraft chooses any buoy whose 
Belative ee (sac |OSe MRO ihe wai Tehrani. lead \ingueuns lem cilia 


craft is then required to fly straight and level at a 
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© 


Y=CPA 


X 


A/C Heading 


CPA = closest point of approach to buoy with 
respect to the A/C heading 


Re. = horizontal ranges to the buoy measured 
ain, © local tons lv2) ancdues 


C = equal distances (assumed) generated by 
flying equal times at a constant true 
airspeed between points 1, 2 and 3 


Note: Buoy coordinates are not known, moreover 
are not required for this problem. 


FIGURE 5 


GROUND SPEED DETERMINATION PROBLEM GEOMETRY 
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constant true air speed for a time itnterVvVa sui) olen uero 
measure from three to eight ranges to Ene (BuO yo in aues om: 
The time between each measurement is short, two to four 
seconds. The time between each range, C, should be exactly 
equal. This Sains be done by electronic time gating. This 
time, C, also represents equal distance between range mea- 
surements. lt is assumed that for short periods of time and 
by electronic time gating, these distances will be equal. 
Any errors are absorbed in the total ground speed error. 


The actual or eames position of the buoy is not 
required as an input to the problem. Also it is not re- 
quired to know the heading or airspeed of the aircraft, 
only that they are constant. The evaluation of C deter- 


mines the ground speed. The use of the Pythagorean theorem 


produces three equations for the three ranges measured; 


2 2 2 
(X+2C) + Y = Ry for range measured at 1 
2 
(X+ C)° + Y* = Ris for range measured at 2 
i 2 2 
X + Y = R3 for range measured at 3 


The magnitude of the ranges only are known and not the direc- 


tion. There are three equations and three unknowns; X, Y 
and C. The values for C are evaluated from the above equa- 
tions for the three ranges; 
2 2 2 
Ra > Zien Ra 
2 


- 


By a similar development C is determined for any--number of 
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ranges measured. Up to eight ranges are used in the de- 
tailed mathematical development included later. 

| Two studies were conducted to determine the feasibility 
of oe ground speed determination method. The first study 
was a sensitivity analysis to determine the effect on ground 
speed error by varying the following; (1) number of ranges 
used, (2) time between range measurement, (3) ranges to the 
buoy and (4) the relative bearing of the buoy. 

The second study was an error analysis on the ground 

speed as affected by errors in the measured ranges. The re- 


sults of these studies are combined in Appendix B. ea 

This method for determining ground speed ts expected to 
have limited application, however, because of the simplicity 
of concept it should be quite reliable to supplement break- 
down of more sophisticated tnertial and dopplar navigational 
systems. Frequently, in fine localtzation, a timed run out 
from some reference to a fix position for attack ts required. 
lt is felt that this ground speed determination method can 
be utilized on the final leg to more accurately determine 
the ground speed and hence a more accurate time to fix on 
that leg. The aircraft will maintain a true course and air- 
speed to that fix position which meets the requirements for 
the solution. 
oe TRUE COURSE DETERMINATION 

lt was felt that the use of the ange measurement to a 


buoy as the only input could provide additional information, 


if applied propery, in evaluating a "drift vector''. Consider 
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im Figure 6-a, an aiteraft flying amb emneoduncmoimcnlc 


North. This true heading is held, at a constant true air 
speed (TAS), for a time ty: Vector Diy represents the true 
air speed times the time flown, ty: The aircraft then 


turns to a true heading of East and holds this heading, at 
the same TAS, for a time te where te = ty: While on each 
of these headings the aircraft determines ground speed by 
methods described in the previous section. A different 
buoy is used for ground speed determination for each of the 
two headings. lf a ad from the Southeast were affecting 
flight, the ground speed for the northerly leg would be 
greater than the TAS and the ground speed for the easterly 
leg would be less than the TAS, [Fig. 6-b]. The triangle 
G-W°D“ shows the second leg geometry without the advance 


Bee? 6—E 
and transfer of the aircraft in the turn. The magnitude 
only and not the Tee of these ground speed vectors is 
known. It should be noted that any initial true heading 
and TAS could be flown and a subsequent right or left turn 
could be executed. Maximum accuracy of the drift vector is 
Be eained when the subsequent heading is perpendicular to the 
initial heading. 

In the Aeon neni don the determination of a drift 
vector, the advance and transfer of the aircraft was elimi- 
nated in further diagrams and does not alter the geometric 
mecuracy of the development. The actual path of the air- 


craft, less -the advance and transfer, [Fig. 7-a], is repre- 


sented by the vectors G, and G, whereas the Vector Diy and 
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nn ne 
e 

Za 

( aE 


/ 
iS Advance and transfer 


of the A/C tn turn 





Dy = Ge - TAS 
= where ty = te Dy = De. 
De = te * TAS 
W, = t - (Windspeed + additional errors) 7 
1 N Wer st 
W. = te * (Windspeed + additional errors) Z 


Ot (Groundspeed for North leg) magnitude only 


cee eee (Groundspeed for East leg) magnitude only 


FIGURE 6 


TRUE COURSE DETERMINATION PROBLEM GEOMETRY 
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(a,b) 


W = drift vector = 


wind vector plus additional 


error 
_ components 


Known Values: 


hyoeassumed position coordinates. (a,b). (c,d) and (est) 
vectors Dy and De, both magnitude and direction 
vectors Gy and Ge, magnitude only 

Unknown Values: 
vectors Gn and G 


Goordinates (x .¥ 


e direction only 
drift vector, W, 


both magnitude and dtrection 


et GUR Go 


TRUE COURSE DETERMINATION SOLUTION GEOMETRY 
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BD represent the ''no wind'' track. The geometry of the pro- 


blem was altered for the purpose of the mathematical solu- 


tion technique. The triangle GeWoD, was Stitt ted Rica =.) 


such that W. coincided with Wy Seine c We = Wie The Wweetor 


DE did not change in orientation thus the relationship be- 


tween the two triangles, GEWoD- and GW. was not changed. 


The resulting form revealed that the vectors Gy and Ge had 


the same terminal point. The known and unknown values in 
the figure are listed at the bottom of Figure 7. 


A pair of quadratic equations were used to solve for 


meme terminal point, (X,Y): 


erie Wenie = G 
(X-c) “+ Gea a Ge 


The magnitude of the drift vector, W, was then solved,» as 
me the determination of the coordinates (X, Y) both the 


coordinates of its origin and terminal point were known. 


wos fx-e) 7 (y-F)? 


The resulting drift vector was then used in conjunction with 


—s 


TAS and true heading to determine an actual course and 


ground speed. 
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V. SONOBUOY LOCATIONS TECHIMO UES 


From the description of the basic sum of squares model 
in section IV, the sum of squares equation for the parti- 
cular geometry in Figure 3 was obtained; 


2 {eae 
{R.-~[(X-A.) +(¥-B.) 9°} 


sg = + 
I 


as 


ie Turther development let 


x. = X~A. Ct) 
Vie = NB. 2) 
then 
1 n 2 2 4 2 
ey rey (3) 


i=] 


Equation (3) was then differentiated with respect to X and 


SSE SSE 
Y and evaluated at Zero, !.€., ae = aY ~ 0 
cee 2a > ee 2s 
aie me ER, (x,+y,) 1k x) (2x) (x. ty.) eee) 
1 Oa a: 7 
= } en) x = (0 (4) 
ee ae ae 
(x, +y,)° 
let (a) be 
R.x. a 


ye 











and let (b) be 


{e435 
x 


| i 
" U=1 : 
Each term of the expression (a) was multiplied, both numera- 


3 1 
tor and denominator, by si 
i 








R ion) R R 
be Lloly 2 
ie n 2 n 2 
ical A Siete i=7 i=1 (1+tan 6.) 
x, ee (1t—))? | 
: ; xe: 
where O. is the relative angle of the buoy at the time the 


range R. was measured, and 


y . 
tan G = e = (5) 


Thus (a) becomes; 


R. , 2 
SS ) Kee OSG. 
ire I 


1 
Nosy (Sec0,) i=1 


N= 


tttias 


The second term, (b) was expanded; 


The two terms (a) and (b), were then combined; 


1 owe ieee 
= ) R.cos0, -X ti ) A, = 0 


i=] i=1 | 


n | n 
X = Ps ) R cos, Saee A. - on © 


nRi— 








A similar development for yy O yielded 
1 ¥ ine 
oe ) R,sin@, - — o) B ; (7) 
i=1 i=] 
The unknowns X and Y could not be solved directly as 0. i's 


Be function of X and Y (see equation 5); therefore, an initial 


value common for all the 0. was estimated, namely Os and 


0 


initial values for X and Y, i.e., X, and 1 Weige determined 
using 8, for each 6, in equations 6 and 7. These values of 
X, and Y, in turn were used to update the values for O., 
namely; 
Ya. = be 
a 0 { eas 
tan — XK, -A, vv Ljee ata 


These values for the 0. were then substituted in the equa- 
mons 6 and 7 to determine X, and Y,. This process was 
repeated until the convergence of X and Y to within accept- 
able values was attained. 

The foregoing iterative process proved to be lengthy. 
The angles 0. were quite insensitive to small changes in 
the positions of X and Y through the iterative process. It 
was felt that an application of Norton's method in two vari- 
ables might provide a more rapid convergence. The method 
involved the expansion of the functions (as determined by 
equations 6 and 7 using @,) in Taylor's Séries in two 
variables about X, and Y,. The expansion would terminate 


Peeccr the linear terms and the remainder tgnored saline 
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The resulting equations would be solved for X and Y and the 
solutions used for an updated expansion. This process would 
be repeated until the values of X and Y were within accept~- 
able limits. 

The sine and cosine of the relative angles 0. as ap- 
pearing in equations 6 and 7 were developed in terms of x. 


and Y; (see equations 1 and 2); 








COSMO. = ae So see 
(sec 6.)° (1+tan 6.)* (ity, /x,)* 
: xen 
- 
1 
(COR 
sin @®. = i SS ere 2 oe eee 


N= 


2 1 2 2 
(csc @.)* (1+cot 6.)° (14x: /y,) 


Yj 


20 24 
(x+y. ) 
The equations 6 and 7 were then cast into the homogeneous 


forms below; 


1 y aa 
315 OS OS {R.———-——- - A.} = 0 _ 
ent (xe tys)# 
ae Yj 
go(X,Y) = Y -— J {R,— > + - B.} = 0 


n ° 
i 


— 
1 Oty )F 7 


The form of the Taylor expansion is 


g(X,Y)=9 (Xo, Volta, (Xo, Vo) (X-Xa) tg, (Xa,¥o) (V-VYo) Remainder 
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Spec | fl Ganllty Ge malic ede. 


2 
| n ae 
gi (X,¥)=Xo-+ J {R, —+ eo ER bs 
| Y= (x t+y.)° i=1 '(x. ty.) f 
+ (4 ¥ ee = 0 (8) 
i=1 (xi ty.) 
and gy, expanded; 
n rae n eave 
1 1 oc 
SPC) Vina ) ames Sa) eile ) Raa ad (XX) 
i=1 (x ty.) ie ERY 
n x? 
+[1-— JY R,—s—s-y HN (Y-Y,) = 0 (9) 
i=1 (xity,) 2 : 


with x; and Ova mUlaluetas a Cine, Yon. 
Equating (8) and (9) to zero yielded the following system 


of linear equations; 


Z 


n ve n x.y; 
Bh Petes | iY 
(x, ty. ee i=] M(x, ty.) 
7 2 
¥ 9 y Sn X 9 y LF 
Sie Noe eye eee aye 
iN eeee M(x. ty,) @ a (xity.) 2 
n oe n 
—_ . —_ ) A, 
n cen 
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2 


i y aia 1 y “| 
[— Ree reese = LU liges R , ——_———_—__] Y 
Y= aos ac i=1 Gee ae 
x — ; 
ell X.Y; el x 
Se tae 1on zo, 
l= calc) eee ay 
n oa n 
- 2 : re y B. 
n n 
: I(x, ye *)3 i=1 
Cramer's rule was applied to solve for X and Y. The coef- 
ficient matrix, H, of X and Y, 
2 
D D oleae " a, aia 
1 3 eee 3 
Die) (xechyey) = So ence cen 
H = = 
2 
n xy n xX, 
Yaa e) k, sr IH; i Ta 7 
i=] Pies ty.) =| (x) +y ) 2 
and the vector of constants, 
, 2 
n ae onl Y. n >a n 
0 1° 1 0 1 1 
D — ) 8. - y ge aaa See } Ce A. 
2 
pT ist (ety, 8 stb OGeyl) Misi (sey) # Mien | 
C = 
<n Xess Yon ae n Y ; n 
0 | 0 | I 
D — ) R. -— } R,———-- +- } — ‘nt 8 
6 ae ee oy oP aaa 
n ay Oc ty ) 72 n a Ce aiy ay) 7o Nn eee cae a n. x 
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were applied in the following manner; 


5 3 
D D Wee fe 
Yow dbo. he) le ee (10) 
Dy D3 Dae b= D,° D, 
DoD 
D D 
] 5 
Pee ee eee i 
; D, D, D, Os D, De 
Do D, 


Equations 10 and 11 replaced equations 6 and 7 in the itera- 
tive process. This method proved highly successful as a 
rapid convergence process. The results of the analysis of 
this process are represented by the figures in Appendix A. 
With the success of the second method an attempt was 
made to interpret the distinctive improvement. The coef- 


ficient matrix, H, was broken down and analyzed; 


2 

a n Y ; il n 7 XY; 

i Z 
1 0 DR see) te x; es 

ia = 

2 

0 1 n x.y n Kr 

le R i’ i iy R 
cette (at yee) iil wx) ty.) 


ee 








Ihe elements of the seconayter newer cme acon cor 


| n R. , n R. 
- } —z——(sin 0.) > —z—5~, (sin@,cosé,) 
a (x. ty.) 2 a (x+y, )2 
I 
n 
R. R 


nN ° 
I . 2 
(simGscose. ) -) ——————_(cos Q.) 
Ve et Con ) 


2 
i=] (x+y. 


The matrix below extracted from the second term of 
the coefficient matrix has a determinant whose value is 


Baiial to zero for any values of O85 


2 
=-sin 0. sinO,cos@, 


2 
sin@.cos COG 0. 


If, in fact, there were no errors between the measured 


* Py 


2 2 
ranges and the actual ranges then ee (x, + y,)= and the 


i 
R. 
factored term, ap would equal one. This implies 
(Cc eyes) 2 
that if there were no error between the measured and actual 


ranges, the second term in the matrix H would go to zero 


leaving; 


lt appears that the existing error between the actual and 


° 


measured ranges causes variations in the value of the second 
term in the coefficient matrix H and therefore provides the 


driving force of the convergence process. 
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VI. GROUND SPEED DETERMINATION 


From the problem description wi nesectronel) set new ciate. 


matical model below for ground speed determination was de- 


veloped. Figure 5 provided the geometry for the problem. 


The particular geometry in the figure indicates three range 


measurements taken at equal intervals. Ground speed, using 


from three to eight ranges measurements was developed in 


this study. The distances C. calculated below are inputs to 


the ground speed determination equations 18 and 19 at the 


end of this-section where i = 3,4,...,8 represents the num- 


ber of measured ranges utilized in the calculation. 


For the three range geometry the known values are Ry» 


oy and R,. It is also assumed that the time between’ range 


foresurements is the same itn each case. The unknown values 


are X, Y and C. The three equations from which C was deter- 


mined are applications of the Pythagorean theorem. 


(X + 20)" + ¥° = RD xX” + ACK + 4c? + ¥? = R® (a) 
2 74 2 2 2 2 2 
(X + Ge en ere eGo eC et) ey) 
é 2 2 2 2 
- X + Y +R, —™ X + Y = (ey) 
3 3 
meetions (a) plus (cc) minus 2(b); s 
2 a 2 2 
eC Kae Cee tsa = ei, 
x” +Y = R, 
-2X° = 4X - 2C° -2Y¥ =-2R, 
Aaes3 
2c" pee acne 
Re aeeR 
a : (12) 


Ly 





The development for four ranges; 


(x + 3C) + Y = R, > X + 6CX + 9C + Y* = R® (a) 
(X20) Sat X20 RE Ge aero re eee ae aL) 
(X + Bo eee Xe eee lean (c) 
x* + Yo = RE xX? # Y= RF (d) 
Equations (a) minus (b) minus (c) plus (d); 
x" + 6CX + 9C° + YY = RK 
-x° = 4cx - Aco - Yo = R, 
-x? = 2cx - c? - y? = R? 
x" +Y = R, 
hee = Ri - Ry RE + Re 
a R, - R, - R, +R, 2 


The development for five through eight ranges was omitted 


but was similar to the two cases above; 





oS (14) 
C, = ; (16) 
c= (17) 





The value C. (S -aeaiustance snmeUm tse ot Nauti Ga lmennolers 


or yards in each of the cases above. 
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The conversion to ground speed was easily accomplished; (1) 


for C. in nautical miles and (2) for C. in yards; 
| 
Cutan) : 
_ i B000(sec),  3o00me, 
(1) GS = THE tsee)~ T(hr) ~~ “Time Knots ve) 


(2) GS = vi tYas? | 3600(sec) 1 (nm) = em ores (19) 
TIME (sec) 1 (Chr) 2000(yds) TIME 





where TIME is the time interval in seconds between each 


meuae measurement. 


— 
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Vit. TRO. COURSE BDETERMINA TRG 


The true course determination used the ground speed as 
calculated in the previous section. Figures 6 and 7 of Sec- 
moon V provided the geometry tor the See e ena Pour e a7 =o 
is annotated for the following mathematical development. 

The coordinates of X and Y were first determined. This was 


done by the simultaneous solution of a pair of quadratic 


equations. 


a) font L a G 
(X-c) + (Y+d)? = Ge 


Expanding these equations; 


a erarea ee by + bo = ef (ay) 

are (C ie divest dic = Ge es) 
Smptract (b) from (a): 

2 (eee A Se Ba (en hes ees es he 


Solve for Y tn terms of X; 


~~ 


2(d-b)Y = 2(a-c)X + G 


y = bare) y-4 ° 
(d-b) 2(d-b) 





Let 
2 2 
Gq? - G + oz + d - a - Be 
pi ese sc F aoe N a 
"we hb oO 2(d=b) 
Then 


-~< 
{! 

~~ 
< 

+ 
” 
oo 
NS 

oO 
See 


Ky 





Substitute the value for Y in equation.i(b) above; 


2 zZ 3 
X selene (R°X+S) - 2d(R°X+S) + d° = Gy 
zz 2 2 2 2 - 
X - 2cX + c + R X + 2RSX + § tee? Gin eee eS 
+d -G. = 6 
. 


Collecting terms; 
2 2 Z 2 2 2 
(14+R )X + (2RS-2c-2dR)X + (c +d +S -2dS-G, ) = 0 


This is a quadratic in X; 


y = 7B +B -HAC on 
2A 
where : 
hoe (ens 


Be 2 hs 2c 2d) 


eee (c'+d°+S°-2dS-G_-) 


Equation (18) is a quadratic form providing two possible 
solutions. The values of the coordinates X and Y determined 
Beeepoth solutions of X provide points which are mirror tmages 
on each side of a line joining the points (a,b) and (c,d) 

in Figure 7-b. The choice of coordinates (X,Y) depends on 
the initial heading of the aircraft and whether a right or 
Bett turn is executed on the second leg. - 


PPoOmethvemeoordinates (xXeyYImethe. dritt vector was then 


determined. The coordinates (e,f) are known. 
2 2 
W= a(X-e) + (Y-f) _ (22) 


The drift vector was then applied to the terminal heading 


Somcetermine the tre eCourse Otetne altecralt fom. aectnarc 


track determination. 
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Vitt. CONCLUSTONS 


The study showed that the method of determining sono- 
BRS location relative to the aircraft with consecutive, 
short-interval range measurement techniques is feasible. 
Moreover the application of the principle of the first 
derivative of the range with respect to time of a buoy near 
the aircraft heading as a means for determining ground speed 
is accurate and relatively insensitive to range errors, num- 
ber of ranges used and the relative angle of the buoy. 

An increase in the number of consecutive ranges mea- 
sured in excess of three, reduces the error in buoy location 
but the reduction is slight. The aircraft should not sac- 
rifice flexibility to acquire additional range measurements 
unless tactics permit. 

In the sensitivity analysis of the sonobuoy location 
technique, the convergence process appeared to break down 
for acute relative buoy angles where the estimate of the 
relative angle lay between the actual relative angle and 
the aircraft heading. Before implementation of this system 
further research should be conducted in this area to pos- 
sibly provide a 360° zone of reliable information on buoy 
location. The model, at present, can provide 360° of re- 
liable information. For acute relative buoy angles, with 


respect to the aircraft heading, however, the estimated 


angle must be greater than the actual angle. 


= 


~ 


lt ts felt that the utilization of buoy ranging in- 


formation for the location of sonobuoys, apart from 


46 





triangulation procedures, has a broad field of application, 
yet uncovered. A preliminary study was conducted to deter- 
mine relative buoy positions independent of aircraft navi- 
gation. Ranges to all buoys in the pattern were measured 


From a single aircraft position. For a pattern consisting 
] 
of N buoys there were oe possible two-buoy combi - 
CZ Niele 
nations. The actual buoy positions were not known thus the 


assumed positions were used. The ranges were plotted in 
reverse from the assumed positions of the buoys generating 


N? 


fe) (N-2)! 
from each two-~buoy combination. From these resulting fixes 


-aircraft fix positions determined by the ranges 


an average position was determined by least squares method. 
An tterative process was then developed whereby the assumed 


buoy positions were adjusted such that the reverse plot- 
Ni 


ting of the ranges would_generate ————————- aircraft fixes 
(2!) (N-2)! 
that were coincident. This did not produce a unique set 


of final assumed positions for the buoys but did establish 
a locus for each buoy in the pattern relative to the air- 
craft _position. {t is felt that further study in this area 
would provide, possibly, a more versatile use of sonobuoy 


— 


ranging. 
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APPENDIX A, SONOBUOY LOCATION ERROR ANALYSIS 


The first study conducted with respect to sonobuoy 
location was to determine the response of the speed of 
convergence of the iterative process as affected by the var- 
iance of; 

(1) The relative angle of the buoy with respect to the 

aircraft heading, 6. 

(2) Varying estimates, 6. of the relative angle. 

‘ol Range trometnemalircrart tomthe buay. 

Figures 8 through 11 show the effect on the number of itera- 
tions required to converge by variance in the estimate of 
the relative angle of the buoy. The actual relative angle 
ranged from 10° to 80° in ten degree increments. The esti- 
mate of the relative angle, in all but one case, ranged 


° to +15° in error of the actual relative angle in 


mmom -—15 
increments of one degree. For the case where the relative 
buoy angle equals 10° the angle estimate, OO» started at 0° 
which was -10° in error of the actual angle. The convergence 


is quite rapid for all estimates, Oo particularly when 


0 6. The convergence process breaks down when the buoy 


— 


is close to the heading of the aircraft and an estimate 06 
between the actual relative angle and the aircraft heading 
is chosen. This can be seen in Figure 8. It should be 
noted that for any relative angle, 8, however small, if an 
eat imate , OO: is chosen which is less acute than 9, the 
process provides rapid convergence. Data for the effect 


on the speed of convergence for variance in buoy range is 
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not included. The effect was felativeol asi nsigme nea ge 
Buoy ranges at 30 nm. and 5 nm. were tested with at most 
Bia iteration difference (in favor of the 30enmasnande min 
the convergence. The data in Figures 8 airtel th Pepire= 
sent a buoy range of approximately 15 nm. and a range error 
Suecero. Lhe convergence process stopped when the buoy 
location error was within .5% of the range to the buoy. 
The second study incorporated a range error to deter- 
mine its effect on final buoy location at the termination 
of the convergence process. Relative buoy angles of 10a 
45° and 80° were investigated. Figures 12, 13 and 14 show 
the total] buoy location error in yards versus the appli- 


cation of the range error in the following manner to the 


three ranges; 


R R R 
no error applied 
+ + + all errors added 
~ - - all errors subtracted 
+ - + other combinations 
- + - other combinations 


There appears to be an approximate one-to-one relationship 
of the average buoy location error to the range error ap- 
plied which is appealing with respect to the small! angles 


utilized between range measurements. 
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lterations for convergence 


[fcereaercens 1ror convergence 
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FIGURE 12 
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Method of application of range error 


FIGURE 13 
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Method of application of range error 


FIGURE 14 


ERROR IN BUOY LOCATION IN YARDS VERSUS METHOD OF APPLICATION 
RELATIVE BUOY ANGLE 80° 7 
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APPENDIX B, GROUND SPEED DETERMINATION ERROR ANALYSIS 


The results of the study conducted on the ground speed 
determination are pictured in Figures 15 through 20. They 
represent the ground speed error in knots versus the range 
error applied, under the varying conditions listed below; 

(1) Number of ranges measured (3,4,...,8). 

(2) Time between range measurements. 

(3 )sesRelatime bDearingmotebuoywto milrcrart headimar: 

[n no case does the ground speed exceed 1.75 knots. 
Moewcase of -the relative buoy angie of zero was not tested. 
[lt was felt that this was a trivial case. The change tn 
horizontal range with respect to time of a buoy dead ahead 
was the ground speed exactly. It was desired to examine 
the accuracy of the method of application of first deriva- 
five of range with respect to time as the relative buoy 


angle increased from zero. 
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FIGURE 15 


GROUND SPEED ERROR {tN KNOTS VERSUS ERROR IN RANGE DUE TO 
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SONOBUOY LOCATION PROGRAM 


START 


INPUT 
Estimate of buoy relative 
eiGiheyc 

A/C coordinates, A.,B. 


Measured ranges, R, 





Calculate initial buoy 
coordinates aot 


O 

fa. ie 
: Oe ~ mie ieee : aie a 

1=1 i=1 

: ie ioe 
ts = =) R. sind — ) B 

i=] i=] 

LOOP 


J iterations 












Calculate buoy coordinates X & Y. 
X and Y = polynomials in terms of 
Ke aniche eG OosGicine hdl ccum tc ative 
convergence process using previous 
values of X and Y to calculate new 
values. 






lteration = 


yes 


OUTPUT 
Buoy locations 
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STOP 
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SONOBUOQY LOCATION SIMULATION PROGRAM 
START 


INPUT 
| A/iGeecoondinetes 24, cb 
2 Actual buoy coordinates, (C,D) 
3 Estimate of relative buoy 


angle, Ors 





Calculate actual 
ranges to buoy 


in = {(C-A,) *+ (SE: 


Establish acceptable error 
test value as .5% of range 


TESTX -005(R,) 
TESTY .005(R.) 



















Calculate initial values of 
buoy coordinates using 0, 


ee fea 
Y = i) R.sin8 = a y B. 
: Nes4 © Mis] 


LOOP 
J iterations 


Calculate buoy coordinates X & Y. 
X and Y = polynomials in terms of 
X and Y. Loop generates iterative 
convergence process using previous 
. values of X and Y to calculate new 
values 
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GROUND SPEED DETERMINATION PROGRAM 
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